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Abstract 

We  prove  central  limit  theorems  (CLT)  for  the  following  problems  in  geometrical 
probability  when  points  are  generated  in  the  [0, 1]'^  cube  according  to  a  Poisson  point 
process  with  parameter  n: 

1.  The  length  of  the  /bth-nearest  graph  Nk^„,  in  which  each  point  is  connected  to  its 
k-th  nearest  neighbor. 

2.  The  length  of  the  Delaunay  triangulation  Z)„  of  the  points. 

3.  The  length  of  the  Voronoi  diagram  Ki  of  the  points. 

We  show  using  the  technique  of  dependency  graphs  of  Baldi  and  Rinot  that  the  de- 
pendence rajige  in  all  these  problems  converges  quickly  to  0  with  high  probability.  In 
addition,  our  approach  leads  to  more  efficient  sequenticd  and  parallel  zJgorithms  for  this 
class  of  problems  on  randomly  distributed  points. 

1      Introduction 

In  a  pioneering  paper  by  Beardwood,  Halton  and  Hammersley  [4]  and  continued  in  Steele 

[11]  it  was  shown  that  the  lengths  L„  of  several  combinatorial  optimization  problems  (the 
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traveling  salesman,  minimum  matching,  minimum  spanning  tree,  minimum  Steiner  tree, 
etc.)  satisfy  laws  of  large  numbers,  when  their  input  consists  of  n  random  iid  points  Xk,  t  = 
1,...  ,n  in  the  cube  [0,1]''. 

It  is  also  believed  (see  for  example  Steele  [12]),  but  yet  unknown,  that  they  satisfy 
central  limit  theorems  (CLT),  since,  while  the  edges  of  the  above  optimal  graphs  are  not 
independent,  the  dependence  "seems  to  be  local".  In  trying  to  make  this  intuitive  idea  pre- 
cise we  succeeded  in  proving  CLTs  for  three  graphs,  which  are  among  the  most  fundamental 
constructions  in  computational  geometry  (see  for  example,  Preparata  and  Shamos  [9]): 

1.  The  length  A'^jt,„  of  the  fcth-nearest  graph,  in  which  each  point  is  connected  to  its  k-th 
nearest  neighbor. 

2.  The  length  V„  of  the  Voronoi  diagram  of  the  points,  which  is  the  collection  of  Voronoi 
polygons  for  each  point.  Given  a  collection  of  points  X ,  a  Voronoi  polygon  around  a 
point  O  is  the  set  of  all  points  in  the  plane  which  are  closer  to  O  than  to  any  other 
point  in  X  (see  figure  1). 

3.  The  length  D„  of  the  Delaunay  triangulation  of  the  points,  which  is  the  graph  defined 
on  the  given  points  with  an  edge  between  them  if  they  are  Voronoi  neighbors  (see 
figure  1).  This  graph  is  indeed  a  triangulation,  i.e.,  a  subdivision  of  the  convex  hull 
of  points  in  triangles. 

Both  the  Voronoi  diagram  and  the  Delaunay  triangulation  are  fundamental  construc- 
tions in  computational  geometry,  since  many  algorithms  for  solving  geometrical  problems 
are  based  on  their  construction.  In  particular,  the  Delaunay  triangulation  is  quite  impor- 
tant for  the  MST,  since  it  contains  the  MST  as  a  subgraph,  i.e.,  points  which  are  neighbors 
in  the  MST  have  to  be  also  Voronoi  neighbors  (see  figure  1).  This  property  leads  to  an 
efficient  algorithm  for  the  MST,  since  one  first  constructs  the  Delaunay  triangulation  in 
O(nlogn)  time  and  then  runs  the  greedy  algorithm  on  its  graph,  leading  to  an  O(nlogn) 
algorithm  for  the  MST  as  well. 


Figure  1:  The  Voronoi  diagram  and  the  Delaunay  triangulation 


We  believe  that  our  results  give  some  partial  insight  on  why  a  CLT  might  hold  for  the 
MST  as  well.  Ramey  [10]  has  attempted  to  prove  a  CLT  for  the  MST,  but  his  approach, 
although  very  interesting,  did  not  succeed  since  he  needed  some  unproven,  but  plausible, 
lemmas  from  continuous  percolation. 


Figure  2:  The  MST  is  a  subgraph  of  the  Delaunay  triangulation 


The  CLT  for  jVi,„  has  already  been  obtained  by  Bickel  and  Breiman  [5],  who  used  very 
complicated  fourth  moment  estimates.  They  wrote:  "Our  proof  is  long.  We  believe  that  this 


is  due  to  the  complexity  of  the  problem.  Nearest  neighbor  distances  are  not  independent". 
We  used  instead  a  simple  geometrical  approach,  which  works  very  well  in  all  these  problems 
and  it  has  the  potential  to  work  for  other  problems  also. 

Our  results  are  established  for  the  Poisson  model,  which  is  asymptotically  equivalent  to 
the  usual  Euclidean  model,  but  more  convenient  to  work  with. 
Poisson  model 

Let  X,,  i  =  1, . .  . ,  Nn  he  the  points  of  a  Poisson  process  with  intensity  n  on  R'',  which  lie 
in  a  d-dimensional  cube,  so  that  Nn  is  a  Poisson  random  variable  with  mean  n. 

The  paper  is  structured  as  follows.  Section  2  contains  the  proof  of  our  main  result  that 
for  all  three  problems: 

limPr{^"-^/^"Ux}  =  .^(x).  (1) 

n-.oo  y/Var[Ln] 

where  <I>(x)  is  the  cumulative  distribution  of  a  standard  A'^(0, 1)  normal.   In  Section  3,  we 
show  that  the  lengths  E[Ln]  of  the  Jt-th  nearest  neighbor  graph, the  Delaunay  triangulation 
and  the  Voronoi  diagram  satisfy: 
E[U] 


lim 


:?  =  Pd.  (2) 

n  J~ 

The  constants  Pj  appearing  in  (2)  have  been  explicitly  computed  by  Miles  [7],  which  is  in 
sharp  contrast  with  the  problems  studied  in  [4],  [11]  with  the  exception  of  the  MST  constant 
that  we  calculated  recently  in  [1].  We  include  for  convenience  simple  direct  computations 
for  the  constants.  Finally  in  the  final  section  we  remark  that  our  approach  leads  to  more 
efficient  sequential  and  parallel  algorithms  for  this  class  of  problems  on  randomly  distributed 
points. 

2      Central  limit  theorems 

In  this  section  we  prove  CLTs  for  the  3  problenns  we  consider.  We  first  outline  the  method- 
ology and  comment  on  its  applicability  to  other  combinatorial  problems. 


2.1     Intuitive  idea 

We  identify  a  certain  geometrical  configuration  An,  such  that  Pr{An}  — ►  1,  whose  oc- 
currence implies  independence  of  the  configuration  at  points  which  are  far  enough.  More 
precisely,  there  is  a  cutoff  distance,  such  that  if  the  event  -4„  happened,  then  determinis- 
tically,  the  configuration  around  a  given  point  is  not  influenced  by  that  of  points  further 
than  the  cutoff  distance. 

For  all  the  three  problems,  the  event  An  is  the  event  that,  in  the  subdivision  of  [0, 1]*^ 
into  p„  =  j-^  equal  subcubes,  each  subcube  contains  at  least  one  and  at  most  e  log  n  of  the 
points  Xk  If  all  the  neighboring  subcubes  around  a  point  are  nonempty,  we  then  show  that 
only  a  finite  number  of  them  determines  the  fc-th  nearest  neighbors  or  the  Voronoi  neighbors 
of  a  point.  Thus,  conditional  on  the  event  An,  the  three  problems  exhibit  "m-dependence" 
(at  the  subdivision  level)  for  some  finite  m.  Applying  the  theory  of  dependency  graphs  from 
[2]  yields  the  CLTs. 

As  mentioned  in  the  introduction,  our  approach  was  inspired  by  the  Ph.D  thesis  of 
Ramey  [10],  who  attempted  a  similar  approach  for  the  minimum  spanning  tree.  In  order  to 
see  how  the  approach  we  described  for  the  three  problems  we  considered  might  generalize 
for  the  MST,  consider  the  event  Bn,k  that  there  exist  2  points  which  are  Voronoi  neighbors 
and  lie  at  some  distance  /  of  each  other,  but  can  also  be  connected  by  a  chain  of  edges 
which  are  all  shorter  than  /  and  the  shortest  such  chain  between  the  two  points  uses  more 
than  it  edges  for  some  large  k.  Such  two  points  will  not  be  neighbors  in  the  MST  and  the 
decision  not  to  connect  them  ts  affected  hy  the  position  of  far  away  points.  Let  B'^^.  be  the 
complement  of  the  event  Bn,k 

Conditioning  in  this  case  on  >i„  fl  5J,  ^  makes  the  problem  m-dependent  with  m  finite 
and  thus  showing  that  Pr{fl„jt}  — *  0,  as  n  — ♦  oo  would  establish  the  CLT  for  the  MST 
(modulo  some  technicalities).  Further  analogies  with  continuous  percolation,  when  d  =  2, 
make  the  above  quite  plausible,  even  though  a  rigorous  proof  has  not  been  submitted. 

To  summarize,  in  the  nearest  neighbor,  the  Voronoi  diagram  and  the  Delaunay  triangu- 
lation,  we  have  established  the  geometrical  configuration  whose  probability  of  occurrence 


tends  to  1  and  which  impHes  finite  range  dependence  and  thus  the  CLT.  In  contrast,  for 
the  MST  we  know  the  geometrical  configuration,  but  we  can  not  show  its  probabiHty  of 
occurrence  tends  to  1.  Finally,  in  the  TSP  and  minimum  matching  problems  we  do  not 
even  have  a  clear  geometric  argument  of  why  dependence  is  local. 

2.2     Dependency  graphs 

The  fundamental  concept  that  captures  the  idea  of  local  dependence  is  that  of  dependency 
graphs.  Introduced  in  Petrovskaya  and  Leontovitch  [8]  and  applied  to  several  problems  by 
Baldi  and  Rinott  [2],  [3]  dependency  graphs  are  defined  as  follows: 

Let  Xa,  a  G  V  be  a  collection  of  random  variables.  The  graph  G  =  {V,  E)  is  said  to  be  a 
dependency  graph  for  Xa,  if  for  any  pair  of  disjoint  sets  A\,  A2  £  V  such  that  no  edge  in  E 
has  one  endpoint  in  Ai  and  the  other  in  A2,  the  sets  of  random  variables  Xa,  a  £  A\  and 
Xa,  a  S  A2  are  independent.  The  following  is  the  fundamental  theorem  that  captures  local 
dependence. 

Theorem  1  ([3])  Let  {Xa,n,  a  €  V„}  be  random  variables  having  a  dependency  graph 
Gn  =  {Vn,En),  n  >  I .  Let  Sn  =  HoeVn  ^<»n,o-^  =  Var[Sn]  <  00.  Let  Dn  denote  the 
maximum  degree  of  Cn  (^^d  suppose  \Xan\  <  Bn  almost  surely.   Then 

\Pr{^ILZll^  <  x}  -  $(x)|  <  32(1  +  V^k'^"':^"^")^.  (3) 


n 


Thus  if  l^"lff^"  —  0  as  n  ^  00 

^^^—^^-^■(0,1).  (4) 

Remark: 

The  above  theorem  quantifies  the  fact  that  enough  independence  (D„  small),  enough  bound- 

ness  (S„  small)  and  enough  variability  (cn  large)  imply  the  CLT. 


2.3     The  CLT 

Let  Lji  be  the  length  of  the  graphs  Nk,n,Dn  and  V„.    We  subdivide  the  cube  [0,  1]''  into 
a  set  C  of  m''  subcubes  with  m''   =    (i^t;^)-     The  subcubes  will  be  denoted  by  indices 
»  =  (»i, .  . . ,  ij),  where  each  ij  takes  values  from  1  to  m. 
We  decompose  the  total  length  L„  as 

Ln  =  E  hn  (5) 

Tec 
where  L'     is  the  sum  of  all  edges  with  both  ends  in  the  subcube  »',  plus  the  sum  of  the 
portion  of  the  edges  within  the  subcube  with  only  one  endpoint  in  the  subcube  i.  Let  A^-' 
be  the  number  of  points  falling  in  the  cube  i  and  An  be  the  event  that  the  subcubes  will  be 
non-empty  and  contain  at  most  e  log  n  points,  i.e., 

>in  =  n^ec{i<^r.n<«'°gM-  (6) 

We  first  show  that 

Lemma  2 

lim  PriAr,}  -  1. 

n— 'OO 

Proof 

The  random  variable  A^^^  is  Poisson  with  parameter  A  =  -^  =  logn.  Thus 

Pr{\  <  ^,^  <  eA}  =  1  -  Pr{Af,„  =  0}  -  Pr{N^,^  >  eA}. 
But  Pr{Nj^  =  0}  =  e-^  and 


<  x_^_^<e-\ 


e-  1        (eA)!  "  (e  -  1)^^  s/X 
where  we  used  Stirling's  bounds  for  the  next  to  last  inequality.  Therefore, 
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-A 


PriK  N^     <eA}  >  l-2e-^  =  1-- 


Since  we  are  using  the  Poisson  model,   the  random  variables  A'''     are  independent  and 

therefore 

Pr{^„}>(l--r'  =  (l--)i^. 
n  n 

Taking  limits  we  obtain 

_      2 

lim  Pr{An}  >    lim  e    ■»«"  =  1 


and  thus  lemma  2  follows.  □ 

We  introduce  now  a  distance  between  subcubes  i  and  j: 

d(ij)  =  maii<r<d{|»r  -  >]}• 

Let 

SlR={UC:   d{lj)<R}, 

i.e.,  S^j^  denotes  the  sphere  of  subcubes  of  radius  R  around  i.  Moreover,  if  A,B  are  two 
sets  of  subcubes,  we  let 

d{A,B}  =  minj^^j^gd(ij). 

The  following  proposition  is  the  heart  of  our  development  of  the  CLT  and  captures  the 
idea  of  "local  dependence"  in  the  three  problems  we  consider. 

Proposition  3  Conditional  on  the  event  A^  having  occurred,  there  exists  a  fixed  number 
R  (depending  only  on  the  dimenston  d)  such  that  if  A,B  is  any  pair  of  sets  of  subcubes 
with  d{A,B)  >  R,  then  the  corresponding  collections  of  random  variables  {L^„,  «  6  A]  and 
{L-r    ,   i  6  B}  are  independent. 

Proof 

Let  X-'  denote  the  set  of  random  variables  {Xy,  v  £i},  i.e  the  collection  of  the  locations 
of  the  points  in  a  subcube.  We  note  first  that  the  point  process  on  [0, 1]''  obtained  from 
the  Poisson  point  process  by  conditioning  on  An,  will  retain  the  property  that  X^  will  be 
independent. 


We  will  prove  the  proposition  for  d  =  2,  where  the  geometry  is  easier  to  visualize,  but 
the  proof  holds  with  no  changes  in  any  d.  Similarly,  for  ease  of  exposition  in  the  k  nearest 
neighbor  graph  we  present  the  case  k  =  I. 

For  the  nearest  neighbor  graph  (with  d  =  2  and  it  =  1),  we  prove  the  proposition 
with  R  =  4.  Note  first  that  if  two  points  Xk  €  i  and  X;  €  j  are  nearest  neighbors,  then 
<^(».  J)  £  2,  since  if  d{i,j)  >  3,  then  a  point  in  a  subcube  at  distance  one  of  i  is  guaranteed 
to  be  closer  than  A'/  (see  also  figure  3). 

Nearest   neighbor 


X.  \     .  X 


Figure  3:  Local  dependence  in  the  nearest  neighbor  graph 


Moreover,  the  decision  which  of  the  points  Xr  belonging  in  subcubes  in  the  sphere  S^^ 
is  the  nearest  neighbor  of  Xk  G  i,  is  in  no  way  affected  by  any  point  Xt  outside  Sj^.  (Note 
that  the  analog  of  this  property  fails  to  hold  for  the  MST,  matching  and  TSP).  But,  given 
the  event  An,  {L^^'  *  ^  ^)  '^  some  function  f{X-^,  i  €  Sa,2}  and  similarly  {Lj^,  i  €  B}  is 
some  function  f{X^,  i  6  Sb.j}-  Since  d{A,B)  >  4  implies  that  Sa.2,  Sb,2  are  disjoint  and 
since  Xj  are  independent,  the  proposition  follows. 

We  now  turn  our  attention  to  the  Voronoi  diagram.  Arguing  as  before,  the  Voronoi 
polygon  of  a  point  Xk  €  »  is  included  in  the  sphere  Sj^j,  since  the  Voronoi  polygon  of  Xk  is 
the  set  of  points  which  are  nearest  to  the  point  Xk  than  any  other  point  Xi.  Thus,  if  Xk  €  » 
and  Xi  6  j  are  Voronoi  neighbors  then  d{i,j)  <  4,  since  the  midpoint  of  Xjt  and  Xi  belongs 
to  the  Voronoi  polygon  of  both.  In  fact,  with  more  effort  we  can  check  that  d{i,j)  <  3,  since 
if  d(i,J)  >  4,  then  any  circle  through  Xk,  Xi  has  to  contain  a  full  subcube  f,  where  a  point 
Xt  lies  (since  all  subcubes  are  non-empty)  and  thus  Xk,  Xi  cannot  be  Voronoi  neighbors. 


The  decision  which  of  the  points  Xr  belonging  in  subcubes  in  the  sphere  S^^  are  Voronoi 
neighbors  is  the  nearest  neighbor  oiXk  €  i,  is  in  no  way  affected  by  any  point  Xr  outside  S^  .. 
Thus,  arguing  as  in  the  nearest  neighbor  case,  /?  =  8  is  sufficient  to  prove  the  proposition. 
D 

Since  we  want  to  apply  theorem  1  we  want  an  upper  bound  on  the  length  of  L-'  ,  which 
is  provided  by  the  following  proposition. 

Proposition  4  Given  the  event  An 

Proof 

Since  every  distance  is  at  most  the  length  ^  of  the  diagonal  of  the  subcube,  we  have  that 

.,n  -       .,n   ^ 

Given  the  event  >4„,  A^;^^  <  e  log  n  and  since  m  =  (j^^)^  we  obtain  (7).  D 

Finally,  in  order  to  apply  theorem  1  we  need  a  lower  bound  on  the  variance  of  L„.  This 
is  provided  in  the  following  proposition. 

Proposition  5   There  exists  a  constant  /j  depending  only  on  the  dimension  such  that 

Var[Ln]>fdn^.  (8) 

Proof 

We  will  use  a  general  technique  from  Ramey  [10],  which  we  illustrate  for  d  =  2.  We 
subdivide  the  cube  [0, 1]^  in  n  squares  of  dimension  -^  with  area  ^  each.  We  identify 
a  particular  configuration,  which  has-a  strictly  positive  probability  of  occurring  and  has 
nonzero  variability.  We  describe  the  construction  for  the  nearest  neighbor  graph.  Let 
e  =  g-T^.  Consider  all  the  subcubes  that  have  the  following  properties: 


1.  They  contain  exactly  2  points  inside  a  circle  with  center  the  center  of  the  subcube 
and  radius  e,  and  no  other  points  in  the  subcube. 
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2.  A  ring  of  points  exist  near  the  outside  of  the  boundary  at  distance  at  most  2e  of 
another.  In  particular,  we  demand  that  all  the  36  smaller  subcubes  of  dimension  i 
that  surround  the  initial  subcube  are  non-empty  (see  figure  4). 


1  -         1 

=      O       I 

1    1  1  1  1  1 

Figure  4:  The  configuration  for  the  nearest  neighbor  graph 


Note  that  no  point  outside  such  a  subcube  has  a  nearest  neighbor  inside  (since  its 
distance  to  a  point  inside  is  at  least  3e  and  there  are  points  outside  within  2e).  Moreover, 
the  two  points  inside  the  circle  are  nearest  neighbors  of  each  other.  Their  distance  d  has  a 
certain  variance  ^  (easily  computable).  Finally,  the  probability  of  a  such  configuration  is 
p  =  e-^"'i^^e-^(^-''')(l  -  e-^'')3«  =  I.SSIO""^  >  0. 

Let  N  be  the  number  of  subcubes  /  =  {i\, . . .  ,1^}  satisfying  the  two  properties  above. 
Then  £'[A'^]  =  np.  Let  F„  be  the  cr— algebra  determined  by  the  random  variables  I'x, . . .  ,z;v 
and  by  the  positions  of  all  the  points  lying  outside  the  N  subcubes.  Then, 

Var[K]  =  Var[E[L^\F^]]  +  E[E[LI\F^]  -  (£;[L„|F„])2]  = 

V^ar[£[L„|F„]]  +  f:[l/ar[L„|F„]]  >  E[V ar^  d,  +  ^  d.|F„]]  = 

16/  xil 

E[Y  Var[d,]]  =  -E[N]  =  crp  >  0, 
i€/ 
where  d,  is  the  sum  of  all  the  nearest  neighbor  distances  of  the  points  inside  a  subcube  ». 

Letting  /2  =  crp  proves  the  proposition  for  the  nearest  neighbor  graph. 
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For  the  t-nearest  graph  the  construction  is  similar  except  we  require  that  there  are  i  +  1 
points  inside  the  circle,  so  that  all  it-nearest  neighbors  are  inside  the  circle  and  its  of  the 
outside  subcubes  contains  k  +  I  points.  . 

For  the  Voronoi  diagram  and  the  Delaunay  triangulation  we  consider  all  subcubes  with 
the  following  properties:  They  contain  exactly  one  point  inside  a  circle  with  center  the 
center  of  the  subcube  and  radius  |  and  exactly  3  points  inside  a  circle  e  and  no  other  points 
in  the  subcube.  The  three  points  in  the  second  circle  are  located  as  follows:  We  subdivide 
the  circle  in  six  equal  circular  regions.  The  3  points  are  located  in  alternate  regions.  As 
in  the  nearest  neighbor  graph  a  ring  of  points  exist  outside  the  boundary  at  distance  at 
most  2e  of  each  other.  The  probability  of  such  configuration  is  pi  >  0.  Note  that  the 
point  inside  has  Voronoi  neighbors  only  the  three  points  in  the  outside  circle.  Fixing  the 
position  of  the  three  points  in  the  outside  circle  there  is  a  certain  variability  ^,  ^  for 
the  Voronoi  polygon  and  its  Delaunay  triangulation  of  the  point  inside,  all  of  which  are 
easily  computable.  Letting  A^i  be  the  number  of  subcubes  /  satisfying  the  two  properties 
above  we  let  Un  be  the  <t— algebra  determined  by  the  random  variables  »i, . . . ,  »Af,  and  by 
the  positions  of  all  the  points  lying  outside  the  N  subcubes  and  inside  the  larger  circle. 
Then, 

Var[U]  =  Var[E[Lr,\Un]]  +  E[E[Ll\Ur,]  -  {E[Ln\Un]?]  = 

Var[E[U\Ur,]]  +  E[Var[U\Un]]  >  E[Var[Yl  d,  +  ^  d.|C/„]]  = 

—E[N]  =  avpx  >  0, 
n 

where  d,-  is  the  sum  of  the  Voronoi  polygon  (Delaunay  triangulation)  of  the  inside  point.  □ 
We  now  have  all  the  ingredients  to  prove  the  CLT. 

Theorem  6    The  length  Ln  of  the  k-nearest  neighbor  graph,  the   Voronoi  diagram  and  the 
Delaunay  tnangulaiton  satisfies  the  CLT 

^lun^Pr{^"-^/^"Ux}.<^(x),  (9) 

n-.oo  yVar[Z,„] 

where  ^(x)  is  the  cumulative  distribution  of  a  standard  normal. 
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Proof 

We  first  establish  that  the  random  variable     ')'•••'■  "-^  given  that  the  event  An  has  occurred 

y/Var[L„] 

is  2isymptotically  normal.  Indeed,  we  define  a  dependency  graph  in  the  sense  of  section  2.2 
between  the  L^^  by  putting  an  edge  between  Lj-^  and  Lj ^  if  d(i,j)  <  R,  where  R  is  defined 
in  proposition  3.  Proposition  3  establishes  the  correction  of  this  construction.  Moreover, 
the  maximal  degree  Dn  of  this  dependency  graph  is  at  most  R  ,  i.e.,  finite  independent  of  n. 
Applying  theorem  1  we  have  that  iKnl  =  m"^  and  D„  <  72'^  from  proposition  3.  In  addition, 

Bn  <  - — ^  °^i' from  proposition  4  and  Var[L„]  >  fdn~^  from  proposition  5. 

Therefore,  since 

\Vn\DlBl        ^i^)[       'y       ]  nl+i 

<   ■ dTTT <  "<i- 


where  uj  is  a  constant  depending  on  the  dimension  only.  Since  m!^  =  j^^  we  obtain 

\Vn\DlBl  Pognrl 

; — ;t  <  ^d i ►  0, 

Var[Ln]'i  ns 

as  n  — »•  oo.   Applying  now  theorem  1,  we  establish  that     ?~   '  "I  given  that  the  event  An 

has  occurred  is  asymptotically  normal  N(0, 1). 

We  now  need  to  show  that  (/„  =     "~    1  "1  is  asymptotically  normal.  But 

Pr{Un  <x}  =  Pr{Un  <  x\An}Pr{An}  +  Pr{Un  <  x\A'JPr{A'n} 

Taking  linnits  and  using  lemma  2  and  the  asymptotic  normality  of  (7„  given  the  event  An 

we  establish  the  theorem.  □ 

Remark: 

We  can  establish  rates  of  convergence  to  the  normal  distribution  by  using  (3).  We  see  that 

^/Var[Ln\  n* 

3      Computations  of  the  expectations 

In  this  section  we  will  find  exactly  £'[L„],  i.e.,  we  will  compute  the  constants  /Jj,  appearing 
in  (2),  using  the  toroidal  rather  than  the  usual  Euclidean  distances,  so  that  there  are  no 
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boundary  effects.  However,  the  boundary  effects  can  be  shown  to  be  negligible  (Jaillet  [6]), 

and  so  the  same  results  holds. 

The  k-th  nearest  neighbor  graph 

Consider  a  given  point  under  the  toroidal  model.  Then  the  length  of  the  nearest  neighbor 

graph  L„  =  "/?„,  where  Rn  is  the  distance  of  a  given  point  to  its  t-th  nearest  neighbor. 

Also, 


E[Rn]=   I      Pr{Rr,  >  r}dr  = 
Jo 

L 


2^ 7i '^'■' 


^•! 


where  cj  is  the  volume  of  the  unit  sphere  in  dimension  d. 
Letting  u  =  n  cj  r''  we  obtain  that 


Letting  n  — ►  oo  we  obtain  that 

E[u\_    1    Aro  +  ;-i) 

--~i;^-^S  0-1)!  <'"> 

The  length  of  the  Voronoi  diagram 

Let  Rn  be  the  length  of  the  perimeter  of  the  Voronoi  polygon  around  a  given  point  O.  Then 

E[Ln\  =  n  ^2"  ■  ^^  ^'''  perform  the  calculation  for  d  =  2.  Let  £'[/»„]  be  the  expected 
contribution  to  Rn  from  a  point  C  of  the  Poisson  process.  Then  £'[/2n]  =  nE[hn]-  Condi- 
tioning on  the  distance  OC  =  r  (see  figure  5),  let  £'[/i„(r)]  be  the  expected  contribution  to 
E[hn]  given  that  OC  =  r.  Then, 

E[hn]  =   /       E[hn{r)]2Krdr. 
Jo 

Moreover,  E[/i„(r)]  =  2/o     Pr{the  length  dl  at  distance  /  from  A  is  part  of  the  Voronoi 

diagram}^/.  But,  Pr{the  length  dl  at  distance  /  from  A  is  part  of  the  Voronoi  diagram}  = 

Pr{no  perpendicular  bisector  cuts  the  segment  OB}  =  Pr{no  point  of  the  Poisson  process 

belongs  in  the  circle  with  center  B  and  radius  BO}  =  e~'**?'  "*"'  ^".  Thus 

72     ,x/2 


E[hn]  =  4^  /       /      re-'((?)'+'')M/dr 
Jo     Jo 
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Computing  the  integral  we  obtain 


which  leads  to 


lim  Ef/inln?  =  4 

n—>oo 


lim 


n-^oo      ,Jn 


2. 


Figure  5:  The  expectation  of  the  length  of  the  Voronoi  diagram 

The  Delaunay  triangulation 

Let  Rn  be  the  sum  of  the  distances  in  the  Delaunay  triangulation  from  a  given  point  O  to 
its  Voronoi  neighbors.  Then,  the  expected  length  of  the  triangulation  is 

Let  h{r)rdr  dO  =  Pr{A  is  a  Voronoi  neighbor  of  O —  .4  is  a  point  in  rdr  d6}.  Then 

r2T     /-x/a 


E[Rn]  =   /       /       rh{r)rdrde. 
Jo     Jo 


In  order  to  compute  h{r)  we  consider  the  circle  passing  through  O,  A  and  a  third  point  of 
the  Poisson  process  B,  which  is  to  the  right  of  the  segment  OA  with  the  property  that  it 
contains  no  other  point  of  the  process.  It  is  well  known  (see  for  example  [9])  that  O,  A  are 
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L(r,t) 


R(r,t) 


Figure  6:  The  expectation  of  the  length  of  the  Delaunay  triangulation 

Voronoi  neighbors  if  and  only  if  no  points  of  the  Poisson  process  lie  in  the  segment  of  the 

circle  to  the  left  of  O^  (see  figure  6). 

Let  Tr  the  radius  of  the  circle  and  L{r,t),   R{r,t)  the  areas  shown  in  figure  6.  Therefore, 

rv/2 


h{r}  =  2  /      e-^<"'''"dPr{T,  <  t} 
Jo 


where  we  multiplied  by  2  because  we  considered  the  point  C  to  the  right  of  OA,  where  we 
could  have  considered  the  case  in  which  C  is  on  the  left.  Moreover, 

Pr{Tr  <0  =  e"^*''"". 

is  the  area  of  the  right  segment  of  the  circle.  U  4>  =  sin'^j^  is  the  angle  ACB,  then  from 

simple  geometry  L{r,t)  =  <i,t^  -  ^^  and  R[r,t)  =  (tt  -  <p)i^  +  ^^.    Computing  the 

integrals  we  find  that 

64 


which  leads  to 


Hm  E[Rn]s/n  =  — , 


,.       E[U]       32 
hm  — -=-  =  — . 

n—oo      y'n  07r 
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4      On  the  relation  of  CLTs  to  efficient  algorithms 

The  approach  we  used  leads  also  to  a  faster  algorithm  to  compute  the  Voronoi  diagram 
and  the  Delanay  triangulation  for  randomly  distributed  points.  From  proposition  3,  in 
order  to  find  the  portion  of  the  Voronoi  diagram  inside  a  subcube,  we  need  only  to  consider 
only  those  subcubes  that  are  at  distance  at  most  R  =  4  from  this  subcube.  Thus,  if  the 
event  An  occurs,  there  are  O(logn)  points  that  contribute  to  the  portion  of  the  Voronoi 
diagram  inside  a  subcube.  Running  an  0{k  log  k)  algorithm  with  it  =  O(logn)  gives  an 
0(log  n  loglog  n)  algorithm  to  find  the  Voronoi  diagram  inside  a  subcube.  Since  there  are 
j^^  subcubes,  we  obtain  an  0(n  log  log  n)  algorithm.  If  the  event  An  does  not  occur,  we 
simply  run  any  O(nlogn)  algorithm.  Since  the  probability  of  An  goes  to  1,  this  algorithm 
runs  in  time  O(nloglogn)  with  probability  1.  In  addition,  since  this  computation  can  be 
done  in  parallel,  the  above  appraoch  leads  to  an  0(log  nlog  log  n)  parallel  algorithm  with 
i-^^  processors. 

In  general,  if  such  an  approach  were  to  work  for  other  combinatorial  problems,  we  would 
obtain  not  only  a  CLT,  but  faster  sequential  as  well  as  parallel  algorithms  for  their  solution. 
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